Remote Synchronization in Complex Networks 
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We show the existence of a novel dynamical state called remote synchronization in general networks of 
coupled oscillators. This state is characterized by the synchronization of pairs of nodes that are not directly 
connected via a physical link or any sequence of synchronized nodes. This phenomenon cannot be observed in 
networks of phase oscillators as its underlying mechanism is the modulation of the amplitude of those interme- 
diary nodes between the remotely synchronized units. Our findings thus show the ubiquity and robustness of 
these states and bridge the gap from their recent observation in simple toy graphs to complex networks. 



Synchronization constitutes one of the most paradigmatic 
examples of emergence of collective behavior in natural, so- 
cial and man-made systems OKI]. Its ubiquity relies on the 
general framework in which it occurs: the interaction between 
two or more nonidentical dynamical units that, as a conse- 
quence, adjust a given property of their motion. As coupling 
between units increases, synchronization shows up as a collec- 
tive state in which the units behave in a coordinate way. Syn- 
chronization phenomena span across many life scales, ranging 
from the development of cognitive tasks in neural systems ||4|] 
to the onset of social consensus in human societies 0]. 

The ubiquity of synchronization in real systems together 
with the recent discovery [EUl] of their real architecture of in- 
teractions has motivated its study when units are embedded in 
a complex network H] . In this way, each unit is represented as 
a node of a network while it only interacts with those adjacent 
units, i.e. those directly coupled via an edge. In the last decade 
many studies have unveiled the impact that diverse interac- 
tion topologies have on the onset of synchrony |[lQr - [l4ll and 
its stability lEMl. In addition, related issues such as that 
of adaptive networks, in which the interaction pattern changes 
according to the degree of synchronization of the system, have 
also attracted the attention of the community Ill9l - l22ll . 

The former studies mainly rely on the study of coupled 
phase oscillators, such as the Kuramoto model I^H], which 
produces globally synchronized systems as a result of the di- 
rect interaction of pairs of adjacent units. However, it has been 
recently found [25] that, for more general oscillator models 
(in which both amplitude and phase are dynamical variables) 
such as the Stuart-Landau model [2], two oscillators, that are 
not directly linked but are both connected to a third unit, can 
become synchronized even if the third oscillator does not syn- 
chronize with them. This novel phenomenon, coined remote 
synchronization, relies on the modulation of the amplitude pa- 
rameter of an intermediary node allowing the passage of infor- 
mation between two of its neighbors for their synchronization, 
even when the former is not synchronized with them. Thus, 
this tunnel-like mechanism is out of reach in ensembles of 
phase oscillators. 



Remote synchronization has been found to occur in very 
specific and simple topologies such as star-like networks in 
which the central node has a natural frequency different from 
that of the leaves. Within this particular setting it was nu- 
merically and experimentally shown [25] that leaves become 
mutually synchronized without the need of the synchroniza- 
tion of the central node. In this Letter, we aim at showing that 
remote synchronization is not limited to the particular config- 
uration of a star-like motif or a tight specification of the node 
frequencies. To this end, we show its ubiquity and robustness 
in the general case of complex networks. 

The dynamical model consists of a network of N coupled 
Stuart-Landau oscillators |2|]. Each node i is characterized 
by two variables, {xi, y%) T , whose dynamical evolution is as 
follows: 
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where a 2 and oji are respectively the amplitude and the (nat- 
ural) frequency of oscillator i when uncoupled. The second 
term on the right accounts for the coupling of the dynam- 
ics of i with its ki first neighbors. The strength of the cou- 
pling is controlled by A (A = in the uncoupled limit) while 
jC = {Lij} represents the Laplacian of the network defined 
as: (i) for i ^ j, = 1 when nodes i and j are connected 
while = otherwise, and (ii) La = —ki. 

To study the synchronization properties of system (Q3 we 
work with the phase variable of each oscillator, defined as 
6i = tan -1 (yi/xi). Then we can measure the degree of syn- 
chronization of any (connected or not) pair of oscillators as: 



(2) 



* mfrasca@diees.unict.it 



where (-} t means an average over a large enough time in- 
terval. We will consider two nodes as synchronized when 
rij > 6, where 6 is a constant threshold that we fix to S = 0.8. 
Nonetheless, we checked that the results presented are robust 
as other values of S yield qualitatively the same outcomes. 

Once two nodes i and j are classified as mutually syn- 
chronized we label their relationship according to the fol- 
lowing three situations: (i) i and j are directly connected 
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FIG. 1. (color online). Evolution of the degree of global synchroniza- 
tion r [panels (a) and (c)] and the number of remotely synchronized 
nodes urs [panels (b) and (d)] for SF (upper panels) and ER (bot- 
tom panels) networks as a function of the coupling strength A and the 
frequency mismatch Acj. In both cases the networks have N — 100 
and k = 2. The other relevant parameters are fixed to a — 1, ui = 1. 
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1), (H) there is a path of mutually synchro- 



nized nodes between them, and ( Hi) none of the former two 
situations hold. While the first two cases are similar, as both 
are examples of synchronization through physical links, the 
third case is analogous to the observed remote synchroniza- 
tion in a star-like network, but in the more general context of 
a complex network. Thus, we define that two nodes i and j 
are remotely synchronized (RS) when they are synchronized 
(Xij > S) and they are not connected by either a direct link or 
a path of synchronized nodes. 

To quantify systematically the extent of remote synchro- 
nization we count the number of RS nodes, defined as the 
number Nrs of nodes that appear RS with at least another 
node in the network. This allows us to introduce the follow- 
ing order parameter: urs = Nrs/N, representing the nor- 
malized number of RS nodes with respect to the total number 
of nodes N. Finally, to quantify the importance that remote 
synchronization has on the dynamics of the system we also 
measure the global level of synchronization through the usual 
Kuramoto-like order parameter: 
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Note that r takes into account the contribution of both syn- 
chronized (fij > S) and not synchronized (r^ < S) nodes. 

As two well-known paradigmatic network topologies we 
have analyzed both Erdos-Renyi (ER) and Scale-free (SF) 
graphs. The former type of networks is characterized by a 
Poisson distribution P(k) for the probability of finding a node 



with k contacts while SF graphs show a power-law distribu- 
tion, P(k) ~ /c -7 . Thus, while in ER graphs most of the 
nodes are close to the mean connectivity (k), SF networks 
display a large heterogeneity in the number of contacts per 
node as revealed from the existence of hubs having k{ ^> (k). 
For their construction we have made use of the model intro- 
duced in i26ll that allows to control the mean connectivity of 
both networks in order to be exactly the same. In the networks 
reported in this Letter the size and mean connectivity are fixed 
to N = 100 and (k) = 2 respectively. Finally, in order to stay 
close to the framework used in [|25ll we have considered a bi- 
modal distribution for the natural frequencies of the oscillators 
so that nodes with high degree (those analogous to the central 
nodes in a star graph) present a larger frequency, ujh, than that, 
uju of less connected (the ones playing the role of leaves in the 
star topology). In particular, we labeled as hubs those nodes 
having ki > k* t2J\\ and assigned them uji = uoh+ii^h while, 
for the rest of nodes uji = uji + . In the former expressions 
£i is a random variable uniformly distributed between -0.025 
and 0.025. 

In Fig. [T] we show the emergence of remote synchroniza- 
tion as a function of the two relevant parameters: the coupling 
strength A and the frequency mismatch of the network hubs 
Auj = Uh — u>i . In particular, we report the behavior of the 
global synchronization, r, [panels l(a)| and |l(c) | and the frac- 
tion of RS nodes, urs, [panels [1 (b)| and |l(d) | for SF (top) 
and ER (bottom) networks. The results are averaged over 50 
different network instances and, for each network we average 
the results over 10 different realizations of the distribution of 
natural frequencies. 

We find that remote synchronization occurs in both types of 
networks in a region of parameters characterized by a strong 
frequency mismatch Auj and moderate coupling A. In fact, 
for low values of the coupling parameter, nodes cannot syn- 
chronize (either in a direct or remote way) as observed from 
the low values of r and urs. On the contrary, for large 
values of A the network is fully synchronized (r ~ 1) and, 
accordingly, urs assumes values close to zero since all the 
nodes are mutually synchronized with their neighbors. As 
panels [T(a)| and [T(c)| reveal, the onset of full synchronization 
requires greater values of the coupling as the frequency mis- 
match increases. In fact, a large frequency mismatch together 
with values of coupling under the threshold for complete syn- 
chronization favors the onset of remote synchronization, as 
observed from the behavior of urs in panels [T(b)| and |l(d)| 

Compared to SF networks, the values of urs in ER net- 
works are greater, thus indicating that remote synchronization 
in ER networks involves both a larger number of nodes. More- 
over, in ER networks the onset of remote synchronization oc- 
curs for lower values of A. ER and SF networks also show 
qualitative differences in the appearance of remote synchro- 
nization: by keeping fixed Auj and increasing the value of A, 
we find that urs in SF networks show two peaks, while for 
ER networks it shows a rise-and-fall behavior. 

In both (ER and SF) cases remote synchronization ap- 
pears as an intermediary state before full synchronization is 
achieved. However, from the analysis of panels [T(a)l and [T(c)| 
one observes that the behavior of r vs. A for a fixed value 



FIG. 2. (color online). Evolution of the average oscillation frequency 
of each oscillator and urs as a function of A for SF (a) and ER (b) 
networks. The mismatch of natural frequencies is Auj = 1.5 while 
the rest of parameters are the same as in FigQ] 



of Auj is qualitatively different in SF and ER networks. In 
particular, in ER networks (panel [T(c)] » a large plateau around 
r ~ 0.5 is set in the region where remote synchronization 
shows up. In this region, the increase of A does not contribute 
to the overall synchronization level, but to a redistribution of 
the average oscillation frequencies of the network nodes. 

This is evident in Fig. [2 where the average values (over 
the simulation time T) of the instantaneous frequency of each 
oscillator are reported along with the parameter urs. The re- 
sults are obtained by increasing A adiabatically from A = 
so that the system starts from a bimodal distribution as dic- 
tated by the configuration for the natural oscillations. As 
A increases, the gap between the two main frequencies val- 
ues of the bimodal distribution decreases until the network 
reaches full synchronization and the nodes oscillate at a com- 
mon frequency. The readjustment of frequencies reveals that, 
for some values of the coupling, the system undergoes a strong 
reorganization, pointed out by the spread of the oscillation fre- 
quencies between the two extreme values. This readjustment 
coincides with the peaks displayed by urs in both SF and ER 
networks. However, the readjustment seems to occur faster in 
SF networks for which the plateau of r is not observed. 

To gain more insight on the relation between the regime 
of remote synchronization and the onset global synchrony we 
now consider the analysis of all synchronized pairs and its par- 
tition into those corresponding to remote synchronization and 
those for which a synchronized physical connection (either a 
direct link or a path of synchronized nodes) exists. To this aim, 
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FIG. 3. (color online). Evolution of the fraction or RS links, r^, in 
SF (a) and ER (b) networks as a function of the coupling strength A. 
The remaining parameters are set as in Fig. Q] 
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FIG. 4. (color online). Evolution of components of physically (a)- 
(b) and remotely (c)-(d) synchronized nodes in an ER network with 
Auj = 2.8, when A is increased by adiabatic continuation from A = 
1.65 (a)-(c) to A = 1.70 (b)-(d). Nodes are colored according to the 
physically synchronized component they belong to when A = 1.65, 
i.e. in (a). The remaining parameters are the same as in Fig[T] 



we define rjij = 1 if nodes i and j are connected either by a 
physical link or by a path of synchronized nodes and rjij = 0, 
otherwise. We then introduce the following parameters: 
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where H(x) is the Heaviside function. Thus, rp and vrs 
represent the fraction of synchronized links due to a physical 
or remote connection, respectively. Obviously, as vp + vrs = 
1, it is enough to report the behavior of trs. 

In Fig. [3] we show the evolution of trs vs. A for several 
values of Auj. The presence of two peaks in the evolution of 
trs in SF networks reveals a similar behavior to that found for 
urs. As Auj increases, the percentage of RS links increases 
and the two peaks shift towards increasing values of A. On 
the other hand, for ER networks the percentage of RS links is 
higher than in SF networks and vrs shows, as in the case of 
tirs, a rise-and-fall trend. The fall in the number of RS links 
points out that the network is able to recruit physical links to 
get synchronized and thus those regions that appeared as RS 
become merged into a single component made of physically 
synchronized links. 

To visualize the progressive substitution of RS links by 
physical ones in the path towards full synchronization we 
show in Fig. [4] for an ER network (with Auj = 2.8) snap- 
shots of both remotely and physically synchronized links for 
two values of the coupling A. In Figs. |4(a)] and |4(c)| we plot 
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FIG. 5. Evolution of components of physically (a),(b), (c) and remotely (d),(e), (f) synchronized nodes in an ER network with Acj = 1.5 and 
(k) = 2, when A is adiabatically increased. Nodes are colored according to the component to which they belong in (a) and, then, represented 
with the same colors in (b)-(f). The networks correspond to the following value of A: (a)-(b) A = 1.75; (c)-(d) A = 1.80; (e)-(f) A = 1.85. 



two networks corresponding to physically and remotely syn- 
chronized links respectively when A = 1.65. In this case the 
network is divided into several clusters of physically synchro- 
nized nodes (the color of the nodes corresponds to the clus- 
ter of physically synchronized links they belong to) and some 
nodes of these clusters appear remotely synchronized with 
nodes belonging to different clusters [as shown in Fig |4(c)| . 
When A is increased to A = 1.70, two of these clusters merge 
together [Fig. |4(b)| through two physically synchronized links 
that connect each cluster to a new node synchronized to each 
of them. Thus, at A = 1.70 two communities, that were re- 
motely synchronized at A = 1.65, fuse into a single one and, 
as a consequence, those RS links between the nodes of the 
two communities reported for A = 1.65 in Fig. |4(c)| disappear 
at A = 1.70 [Fig. |4(d)) . A further example of the merging 
of RS clusters is shown in Fig. 0. We consider an ER net- 
work with Auj = 1.5 and < k n >= 2, when A is increased 
with continuation from A = 1.75 to A = 1.85. For A = 1.75 
the network is divided into four main components of physi- 
cally synchronized nodes plus some small communities and 
isolated nodes ( Fig. | 5(a)| ). The analysis of the components of 
RS nodes (Fig. |5(d) ) reveals that the there are RS links be- 
tween the four communities. In fact, increasing the coupling 
to A = 1.80 three of these communities merge (Fig. |5(b)| and, 
correspondingly, the RS links between these communities dis- 
appear (Fig. |5(e)| ). Finally, a further increase ofA(A = 1.85 
in Fig. |5(c)| ) leads to the aggregation of the fourth community 
(the bigger one) to the previous ones. Also in this case, amost 
all the RS links disappear (Fig. |5(f)| ) and very few RS links are 
observed for A = 1.85. 

Summarizing, in this Letter we have shown that remote syn- 
chronization can be observed in general complex networks. 
This phenomenon relies on the mutual synchronization of 
pairs of uncoupled nodes. Each remotely synchronized pair 
of nodes are thus physically connected through an intermedi- 
ary node (not synchronized with them) or a sequence (path) of 
intermediary nodes. Although its original formulation was re- 



stricted to a rather particular and simplistic setup, a star graph, 
here we have extended its formulation to the general case of 
complex networks such as Erdos-Renyi and Scale-free graphs 
of coupled oscillators having amplitude and phase as dynam- 
ical variables. The addition of amplitude as a dynamical vari- 
able, in contrast with the typical framework of networks of 
coupled phase-oscillators, provides the observation of remote 
synchronization, a mechanism anticipating synchronization 
by physical links in networks with heterogeneous distribution 
of natural frequencies. 

It is important to note that remote synchronization is not 
observed for ensembles of phase oscillators as its underlying 
mechanism is the modulation of the amplitude of those inter- 
mediary nodes allowing the information transfer between un- 
coupled pairs of nodes. On the contrary, when similar settings 
are applied to phase oscillators a different phenomenon is ob- 
served, namely that of explosive synchronization lfl4ll in which 
the typical second-order synchronization transition transforms 
into a first-order one. In its turn, remote synchronization ap- 
pears as a rather robust state prior the onset of global syn- 
chronization since for a wide range of coupling strengths al- 
most all the nodes are remotely synchronized with, at least, 
another one while the level of global synchronization remains 
small. Thus, our results open the door for experimental obser- 
vations of this novel state in which the existence of a synchro- 
nized pair cannot be associated to a given physical interaction 
through a single link of the network and highlight the impor- 
tant difference between the real (i.e. associated to physical 
links) and the functional (i.e. emerging from synchronization) 
connectivity of a network. 
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